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Abstract 

In this paper we study the density noise correlations of the two component Fermi gas in optical 
lattices. Three different type of phases, the BCS-state (Bardeen, Cooper, and Schieffer), the FFLO- 
state (Fulde, Ferrel, Larkin, and Ovchinnikov) , and BP (breach pair) state, are considered. We 
show how these states differ in their noise correlations. The noise correlations are calculated not 
only at zero temperature, but also at non-zero temperatures paying particular attention to how 
much the finite temperature effects might complicate the detection of different phases. Since one- 
dimensional systems have been shown to be very promising candidates to observe FFLO states, we 
apply our results also to the computation of correlation signals in a one-dimensional lattice. We 
find that the density noise correlations reveal important information about the structure of the 
underlying order parameter as well as about the quasiparticle dispersions. 
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I. INTRODUCTION 



Recent studies in the experiments on ultracold Fermi gases have shown great potential 
for elucidating long-standing problems in many different fields ofphysics related to strongly 
correlated fermions. For instance, in recent experiments [l], y, 0, 0, 0] spin-density im- 
balanced, or polarized, Fermi gases were considered. Among other things, such systems 
make it possible to study pairing with mismatched Fermi surfacesjpotentially leadin g, to 
non-BCS type pairing such as that appearing in FFLO-states 0, 3 or BP-states [ol, lol 
(Sarma-states). These possibihties have been considered extensively in condensed-matter, 
nuclear, and high-energy physics 

In optical lattices it is possible to study many different physical problems with close 
analogs in the field of solid state physics. However, in contrast to many solid state sys- 
tems, ultra-cold gases in optical lattices provide a very clean environment. In other words, 
these systems have very few imperfections and if imperfections are of interest, they can be 
introduced in a more controlled manner. In addition, since optical lattices are made with 
lasers, the lattice geometry is easy to vary [l2|, H, l3, lH| by changing the properties of the 
intersecting laser beams. For these reasons optical lattices enable one to investigate various 
quantum many-body physics problems, such as those related to Mott insulators, phase co- 
herence, and superfluidity. Indeed, the possibility of a superfiuid alkali atom Fernai gas in 
an optical lattice has been recently studied both theoretically [l^, 17, [3, 19, 20, [21I, 22 



as well as experimentally [23|. Furthermore, Feshbach resonances, molecules of fermionic 
atoms, and p-wave interactions in a lower dimensional fermionic systems have been studied 
using optical lattices [2I, 25, 26l . 



Density-density correlations tell us how strongly the atomic densities at different positions 
are correlated. Often it is useful to focus on the deviation of the atomic density from its mean 
value and subtract the background away. In such cases one studies noise correlations i.e how 
the density fluctuations at different positions are correlated. Measuring noise correlations 
is a promising way to observe different phases in an optical lattice since while densities can 
be very similar for different phases, the noise correlations can still be very different. As 
an example, one can mention the famous phase transition between the superfiuid Bose gas 



and the Mott insulator |27l . |28| ]. In the Mott insulator phase one can see the Bragg peaks 
in the noise correlations, but the noise correlations vanish for a Bose-Einstein condensate. 
Experimentally noise correlations have been used as an indicator for different phases in 



optical lattices for bosons [2^] and to observe fermionic anti-bunching for ideal fermions [29 
Also, noise correlations have been used to detect pairing correlations in an interacting Fermi 
gas in a harmonic trap (30| . 

In this paper we study noise correlations in a cloud of ultra-cold two component Fermi 
gas at finite temperatures. We are motivated by the fact that via noise correlations one can 
see subtle correlation effects which are not visible in the lower order correlations functions. 
Not only are the noise correlation signals very different for a superfiuid and a normal Fermi 
gas, but noise correlations also differ between different types of paired states. While we 
compute density-density correlations also at zero temperature, special attention is payed to 
finite temperature calculations, since finite temperature effects have not been extensively 
discussed before. We find that although finite temperature effects do smooth out some 
sharp features present in zero temperature calculations, qualitative and clear quantitative 
differences between phases can still exist at higher temperatures. Also, we show how regions 
of gapless excitations are reflected in the density-density correlations. For instance, gapless 
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excitations show up as absence of correlations for certain momenta. Moreover, they are also 
reflected in the value of the maximum correlation peak height, this quantitative signature 
persists at finite temperature as well. 

Reducing the dimensionality of the system from three dimensions into one has been shown 



to favor the FFLO-type modulated order parameters in free space [31|, l32|, |33| as well as in 



a lattice [3J]. This effect of dimensionality, together with the general tendency of a lattice 
geometry to favor FFLO state due to nesting of the Fermi surfaces ^2|, makes one dimen- 
sional optical lattices a promising system to study the FFLO state. For this reason, we 
apply our general results on density-density correlations also to the one-dimensional (ID) 
system. Very recently a similar one-dimensional problem was also discussed at zero temper- 
ature using the DMRG algorithm [35|. Here we compute the density- density correlations 
of the one-dimensional system at representative points of the finite temperature phase dia- 
gram of a polarized system. The noise correlations in a one-dimensional system turn out to 
show differences between different states which are straightforward to interpret and contain 
useful information on the underlying pair-wavefunction as well as on the structure of the 
quasiparticle dispersions. 

This paper is organized as follows. In Sec. [TTl we discuss the physical system and present 
the Hamiltonian of the system. In this section different paired phases are also discussed In 
Sec. Hn] we proceed to compute the noise correlation functions for different paired phases 
at zero temperature and in Sec. [IV] computations are generalized to non-zero temperatures. 
In Sec. |V] noise correlations in a one dimensional lattice are discussed while making a clear 
connection with the computed ID phase-diagram. We end with some concluding remarks 
in Sec. ED 




II. THE HAMILTONIAN OF THE SYSTEM 

We assume in our calculation a three dimensional cubic optical lattice. The system is 
composed of a two component Fermi gas where components are different hyperfine states (for 
concreteness we assume ^Li atoms). In terms of the field operators \l/(T(r) the Hamiltonian 
of the system is given by 

drrfr'^{(r)^|(r')(7(r,r')^i(r')^T(r)- J2 

E / d^^iir) + ^.(r)+ 

drdr'^{r)^{r')g{r,r')^^{r')^^{r), 

where h = /i/(27r) and h is the Planck constant, fi^ is the chemical potential of the compo- 
nent a and the lattice potential is given by V^(r) = sEr Yl^=i sin^{kxj), where s is the lattice 
depth and Er = h'^k'^/{2m) is the recoil energy {k = Tc/d and d is the lattice constant). In 
the usual way the interaction between atoms is modeled by a contact potential 

g{r,r) = 5(r-r), 

m 



(1) 
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where a is the s-wave scattering length. The number operator of the component a 



is expressed in terms of the density operator n^{r) = ^J.(r)^o-(r). In our case the field- 
operators are fermionic which implies fermionic equal time anticommutation relations 

{^'t (r), ^^(r')} = ^l(r)^^(r') + ^^ir')¥^ir) = 5,^5(r - r') 
{^l(r),^^(r')} = {^a(r),^,(r')} = 0, 

where Sap is the Kronecker delta, and 6{r — r') is the Dirac delta-function. 
A. The Hubbard model 

In a sufficiently deep lattice we can expand the field operators in terms of the well localized 
(lowest band) Wannier functions as 

"^^rir) = ^Wi{r)c^,i, (3) 

i 

where ifj(r) is Wannier function centered around a lattice point i, and Ccr,i is a fermionic 
annihilation operator which annihilate fermions of component a at site i. 

In assuming that only the lowest band states are occupied we are assuming temperatures 
which are much lower than the bandgap. We can estimate the energy of the vibrational 
levels by approximating lattice wells with harmonic oscillators and in this way find that 
temperatures should satisfy a criterion 



, ^ hn 2sEr 

where ks is the Boltzmann constant. If, for concreteness, we assume ^Li atoms in a cubic 
lattice with a lattice spacing 505 nm, the above condition implies temperatures much below 
5/iK for a lattice depth of lOEr. 

It is assumed that only tunneling between the nearest neighbors is of importance. This 
assumption of tight binding is reasonable when the Wannier functions are well localized. 
In other words, overlap integrals between the next nearest neighbor Wannier functions are 
small compared to overlap integrals between the nearest neighbor Wannier functions. This 
always happens when the lattice is deep enough 3^. Using this approximation one finds 
the Hubbard Hamiltonian [3j 



n n 

\ (5) 

{n,m)x (™,'Tt>y (n,m)z ) 
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where {m,n)a means a sum over the nearest neighbors in the a-direction and hopping 
strength is defined by 

Ja = - drw*^{r) I — ^— + sE^^sin^(fcxi) j Wn±d£c(r)- 

Furthermore, the dominant on site interaction couphng strength is given by 

U = ^ [ rfrKWr. 
m J 

In our calculations we have approximated the Wannier functions with a harmonic oscillator 
ground state in each potential well. 

The most stringent criteria for temperatures follows from the fact that critical tempera- 
ture for the superfiuidity is typically much less than the Fermi energy, which in turn is less 
than the bandwidth. This means that in the regime of greatest interest, the temperature 
should be of same order or less than the tunneling energy scale J. 



B. The ground state ansatz at zero temperature 

On a mean-field level in the homogeneous system the ground state ansatz at zero temper- 
ature which includes the possibility of the BCS-state [38)], the breach pair(BP)-state jo]. 



and the plane wave FFLO-state [7|, l8j can be expressed as [39 
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I'^GS) = n 4,-k+q n 4,k n (^^.q + ^k,q4,kCi-k+q)|0), (6) 
keGs keG2 keGi 

where |0) is the vacuum state. In area Gi both quasiparticle dispersions i^f^k.q, -E'j„k,q are 
positive, in area G2 -E|,k,q > and -Et,k,q < 0, and in area G3 -Et,k,q ^ and -Ej.k.q < 0. 
Area Gi is populated by pairs with a quasi-momentum q. Area G2 is populated by atoms 
of the component t in the quasi-momentum states k G 6*2 while the area G^ is populated 
by J,-atoms in the quasi-momentum states — k + q. By diagonalizing the mean-field theory 
one finds the quasiparticle dispersions 



Si/,,.,, = \l I attf^) % T 2±^|i^. (7) 

where A is the pairing gap and the single particle dispersions are 

3 



e<T,k = ^ 2 Jq(1 - cos{kad)) - /i^. 



a=l 



The quasiparticle amplitudes Mk,q and fk,q are given by 

et.k + Ci-kH 




V[(eT,k + ei,-k+q)/2]2 + A2 / (8) 
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FIG. 1: [Color online] The phase diagram of Fermi gas in a 3D-lattice 22] ■ Colors (or shading): 
BCS/BP states=blue or dark gray, FFLO=yellow or light gray, phase separation=red or gray, and 
normal gas=white. The parameters where such that the average filling fraction (n-f + n^)/2 = 0.2, 
J = {)mEr, and U/{2J) = -1.86. 



Based on the properties of this ansatz, it can be classified further. When the wavevector 
q 7^ 0, the ground state is called the FFLO-state which breaks the translational symmetry. If 
q = and both quasiparticle dispersions -E'|,k, -f'i.k are positive for all k-vectors, the ground 
state is called BCS-state, and finally if q = and E^^y^ < for some k-vectors, the ground 
state is called the BP-state. This last case involves phase separation in the momentum 
space. 

In Fig. [1] we show on example phase diagram of the polarized system in a three- 
dimensional lattice [22|. At low temperatures, phase separation is favored for low polar- 
ization, while for larger polarization FFLO state is energetically favorable in substantially 
large parts of the phasediagram. At non-zero temperatures BP-state can be energetically 
stable and occupy regions of the phasediagram close to the critical temperature. Although 
we do not show it in the figure, for low polarization quasi-particle dispersions of the BP- 
state can be gapped while for higher polarization they are gapless. Very close to the critical 
temperature the dispersions are always gapless for the polarized system. Note that only BCS- 
normal phase separation was considered, not BCS-FFLO coexistence which could make the 
parameter window where FFLO occurs even wider. 

III. NOISE CORRELATIONS 



Our aim is to compute the noise correlation functions which are defined by [2 



G 



(n„(r,t)n/3(r',t)) - (n„(r, t)) (n^(r', t)) 



(^t (r, t)i>Ur\ t)i>p{v\ t)^,(r, t)) + (5„^5(r - r') (ri„(r, t)) - (fi„(r, t)) (n^(r', t)). 



(9) 
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where a and (3 are component indeces, r and r' are the positions, t is time. The term with 
Dirac's delta function is due to the normal ordering of the density-density expectation value. 

Often in experiments a gas which has been trapped is released and the gas expands. When 
collisions during the expansion can be ignored and the expansion is ballistic, the density- 
density correlations after the expansion reflect correlations in momentum space at t = 0. 
Therefore, correlations can be computed using the wave- function in momentum space prior 
to expansion. Positions in real space after time t are related to k- vectors trough r = tk/m. 
In other words at long times (the saddle point approximation of the free evolution amounts 
to the limit t ^ md? /Ti) 

^,(r,t)~e-^«W-^% 
where Rj = d{ix,iy,iz) {in are integers) is a lattice vector and 

Thus the densities become 

n.(r,t) = (n„(r,t)) = ^(t)^ ^ e^^W'^- (cl/«,), (10) 

where Rjj = Rj — Rj and A{t) is a time dependent scale factor which depends on the 
Wannier functions. The correlation functions in Eq. ([9]) now become 

i,j,m,n (11) 

+ SapSir - r'){ha{r,t)) - (^^(r, t))(fi/3(r', t)). 

A. The noise correlations of the BCS-state at zero temperature 

We review here the noise correlations of the BCS-state. The usual unpolarized BCS 
ground state is extracted from the ansatz in Eq. by setting both chemical potentials to 
be equal and q = 0. For the BCS-state both quasiparticle dispersions are positive. Because 
the BCS-ground state is presented in momentum space, for simplicity we express the densities 
in Eq. flTU]) in terms of the operators in the momentum space 

(n.(r,t)) = ^ ^ e^«W-^»^e-^'^^-^'+^^-^^(cl,,^c„,k.) (12) 

i,j,ki,k2 



and similarly for the noise correlations in Eq. (ITT 



i,j,m,n 



ki,k2,k3,k4 



giQ(r)-Rij+tQ(r')-R,„„g-tki-Ri+ik2-Rj-jk3-R,„+ik4-R„ 



^a,ki^/3,k3^/3.k4Ca,k2) 1 



+ 6ai3S{r - r'){ha{r,t)) - (^^(r, t)) (^/^(r', t)). 
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where M is the number of lattice sites. Now one needs the expectation values 
(cL,ki4,k3^/3,k4Ca,k2) and (cj^^ki^a.ka)- For the latter term one gets 

(^BCs|cI„kiCa,k2l*iJC5) = ^kikal^kiP (14) 

and using this the densities in Eq. f|T2|) can be expressed as 



R(r,t)) = ^5:e^(«W-'^)-^^^>,r. (15) 

Now the number of particles is the same for both the components as it should be. When 
the lattice is large, i.e. when the lattice size M ^ 1, the exponential terms only add 
constructively when (^(r) — k = and we can approximate the densities as 

(ri„(r, t)) = Aitf 5^ 5 - ^kj 1^;^^. (16) 



Non-vanishing four operator expectation values are 

(^BC5|c|,kcI -k'^i -k'Ct.kl^BCs) = 4k'|^^kP + (1 " 4k' ) I ^k' H ^^k P 
(*BC5|c}{i),k4a),k'^T{i),k'CT{i),k|*i3Cs) = (1 - 4k' ) I ^^k' H ^^k P 

and combining these results and taking the limit of large lattices, we find 
Gn(r, r') = G,,(r, r') = A{t) In^l^kr^ ( r - — ) 5 f r' + 



(17) 



a.(r,r') = -A(t)^J]|t;k|'5 r- 



k 



m 




m 

htk\ . / , mk'\ (18) 



+ (5(r-r')(n,(r,t)), 

where k = k + Yla=i '^naT^Xa/d and k' = k + X]a=i ^rriaTcXa/d, where k is a lattice momen- 
tum, Xa are orthogonal unit vectors and and rrii are integers. Because l^kPl'^kP ~ 
the noise correlation between different components vanishes with vanishing gap A and 
is thus identically zero in the normal state. One can see that the noise correlations in the 
BCS-state differ crucially from this. Because in the BCS-state momenta k and — k are 
correlated, the noise correlation signal is pronounced only if 



+ r' = E 



^ h2ni7TtXi 



i=l 



md 



The continuum result [27|] is similar, but peaks appear when r + r' = 0. In the continuum 
limit our result reduces to the known result. 

On the other hand correlations of the single component G^^ and G^ show a hole when 



, 2hniT(tXi 
^ md 
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FIG. 2: [Color online] In (a) we show a cut of the BCS noise correlations between components in 
the z = plane. In (b) we show the integrated BCS correlations between the components. In (c) 
we demonstrate the noise correlation of a single component in the plane z = 0. As parameters we 
used the average filling fraction (n^ -|- ni)/2 = 0.30, polarization P = (n^ — ri|)/(n| -|- nj^) = 0.0, 
A/(2J) = 1.03 and U/{2J) = —3.0. All these examples are calculated at T = 0. In (a) and (b) we 
choose r' = — r and in (c) we choose y' = y, z = z' = 0, and x' = x + ht27r/{md). Color-coding is 
such that light colors imply high peaks and dark colors low. 



and rii for at least one i [40|]. The reason for this is the fermionic antibunching related to 



Pauli blocking i.e. two or more identically fermions cannot be in the same momentum-state. 
This fermionic antibunching was also measured in the noise correlations of an ideal Fermi 



29 



gas 

In Fig. [2] we demonstrate the BCS noise correlations at T = 0. In Fig. [2] (a) we show a 
cut in the z = plane of the BCS noise correlation between the components while in Fig. [2] 
(b) we plot the column integrated BCS-state noise correlations 

Gi^l{x,y) = J dz Gi^^i {x,y, z, —X + htqx/m, —y + htqy/m, —z + htqz/m) . 

As one can see from the Fig. [21 the noise correlations reflect underlying lattice structure 
and strongly non-spherical structure of the Fermi surface. We can see also from Fig. [2] (a) 
and (b) that the peaks are highest near the Fermi surfaces where e-f,k + e|,_k = 0. In the 
weakly interacting BCS limit correlations are strongly peaked at the Fermi surface, but the 
distribution of the peak heights becomes broad in the BEC limit, and the signal approaches 
the interference signal of the Bose-Einstein condensate. 

In Fig. [21(c) we show the BCS state correlations of a single component in the z = plane. 
As is clear, the BCS-state noise correlation of a single component shows similar antibunching 
behavior as the ideal Fermi gas noise correlation. If r = r' -|- Yl^=i '^hriiTxtXi/ [md) and and 
nj 7^ for at least one i, then the noise correlation of a single component shows the holes. 
However, the BCS-state noise correlation of a single component does also show the bunching 
peak when r = r'. This result differs from the ideal gas result where such a peak is absent. 



B. The noise correlations of the BP-state at zero temperature 

The ground state in Eq. ^ is the BP-state when q = and one of the quasiparticle 
dispersions ^^^^k is negative for some values of the momentum k. This implies that in the 
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BP-state \Sfi\ = I /if — /i J > 2 A. The BP-state is a gapless state, because one of the 
quasiparticle dispersions changes its sign when the momentum varies. The BP state at 
zero-temperature is expressed as 



I'^BP) = n n 4,k n K+t^ka|,kci-k)io)- 

kGGs kgG2 keGi 



(19) 



In the area Gi both quasiparticle dispersions -Et.k.q; -^i,k,q > 0, in the area G2 -E'l.k > and 
-£'T,k < 0, and in the area G3 i?|_k > and ^ < 0. 

Now we can calculate the required expectation values similarly as before. We find that 
the densities in the BP-state (in a large lattice) are 



[n 



t)) = A{tr ^(i^T(i),k) + 0{E,,^)e{Eu^.) W\ 



(20) 



and non-vanishing expectation values in the noise correlations in Eq. (fT3|l are in turn given 
by 

(^Bp|c|,kCl,-k'Ci,-k'CT,k|^Bp) = 

(1 - 5kk')[(l - ^(^T,k))(l - ^(^i,k')) + ^^(^T,k)^^(^i,k)^(^T.k')^^(^i,k')l^k|>k'r 
+ ^(ET,k)^(i5i,k)(l - ^^(i5^i,k'))l^k|' + ^(i?T,k')^(^i,k')(l - ^(^;T,k))|^^k'|'] 

+5kk'^(^T,k)^^(^i,k)it^kr (21) 

(^Bp|c|{i),k4{i),q^T{i),k'CTa),k|^Bp) = 

(1 - 5kk')[(i - ^^(^Ta),k))(i - 0{E^a),^,)) + ^(E^,k)^(i5i,k)^(i5^T,q)^(^i,k')l^k|>kf 
+ eiE^,^)eiE^,^){i - 9{E^a),^,))\v^\' + eiE^y)eiE^,^,){i - eiE^,a)^))M']]- 

Subtracting the densities cancels out most of these terms and again for large lattices we get 
strongly peaked behaviour and 



6 r- 



m 

ntk 

m 



m 



6 r'- 



htk' 



m 



e(i?T,k)^(i?i,k)|^k| 



(22) 



+ (l-^(E,,k)) 



+ 5(v-v'){n,(v,t)). 



In the correlation between different components there are now areas where peaks are absent. 
The reason for this is that the BP-state has phase separation in the momentum space. 
Because in the BP-state there is no pairing between momenta k and — k, if £"1 k < or 
-E;,k < 0, in contrast to the BCS-state, there are no correlations between the points 



, \ ^ h2niTitXi 

i=l 



rad 



in this region. 
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FIG. 3: [Color online] In this figure we demonstrate examples of the BP noise correlations between 
the components at T = 0. In (a) we demonstrate a cut in the z = plane while in (b) we 
show the column integrated BP noise correlation. Parameters we used were A/(2J) = 0.318, 
(/U^ + ^|)/2.0 = 2.18, and 5fi = fi-^ — fi^ = 1.14. In figure (a) the correlations vanish in the dotted 
area. We have chosen r' = — r and color-coding is such that light colors imply high peaks and dark 
colors low. 



In Fig. [3] we demonstrate the BP-state noise correlations at T = 0. In Fig. [3] (a) we 
have plotted the BP-state noise correlation at the z = plane while in Fig. [3] (b) we show 
the corresponding column integrated noise correlation. From Fig. [3] (a) one can see clear 
areas without correlation peaks in the BP noise correlation due to the fact that these areas 
are populated only by atoms of the majority component which are unpaired. From the 
integrated correlations in Fig. [3] (b) we can see that the peaks are now relatively low in the 
area between the Fermi surfaces. This differs substantially from the integrated BCS noise 
correlation (see Fig. [2] (b)). 

Results in this section are intended mainly as a reference since the BP-state is not ex- 
pected to be the ground state at zero temperature. Depending on polarization the ground 
state at zero temperature can be the BCS-state, phase separation [4l|, or the FFLO type 
state |22| . However, at non-zero temperatures the BP-state can be the minimum of the free 
energy [42|. 



C. The noise correlations of the plane wave FFLO-state at zero temperature 

The FFLO ground state is presented in Eq. ([6]). Contrary to the BCS-state and the BP- 
state, in the FFLO-state the parameter q 7^ and the translational invariance is broken. In 
this case the order parameter is given by A(Rj) = Aq exp(iq-Rj), where Rj is a lattice vector. 
In the FFLO-state the other Fermi surface has been effectively shifted by the wavevector q. 
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The densities in the FFLO-state are (in the hmit of large lattices) given by 
(n,(r, t)) = A{tr E ^ - ^) - ^(^T,k,q) + ^(^T,k M'] 



(23) 



In the same way that we computed the noise correlations in the BCS- and BP-states, we 
find the noise correlation in the FFLO-state as 



G,,{v,v')^A{trJ24'^ + 

k 

G,,{r,r') = -A{trJ2 



ntk 



m 



m 



ht{k - q) 



5 r' 



htk' 



6\t- 



m 

ntk 

m 



m 



^(^T,k,q)^(^i,k,q)|^ik,q|>k,q| 



^(£^T,k,q)^(£^i,k,q)kk,q|>k,q| 



Sir'- 



htk' 



m 



(24) 



+ (l-^(Et,k,q)) 



+ (5(r-r')(n|(r,t)) 



G,,{r,r') = -A{trJ2 



6\r- 



ht{k - q) 



m 



m 



^(Et,k,q)^(i?i,k,q)|^k,qr+(l-^(i5i,k,q)) 



+ 6{r-r'){h^{r,t)). 



Because in the FFLO-state momenta k and — k -|- q are correlated, the points 

3 

r 



, Ut h2nintxi 
' + r' -\ q = > 

m '-^ 

1=1 



are now correlated in the expanded cloud. Since the FFLO-state arises only when there is 
some polarization in the system, at T = the FFLO-state is always a gapless state i.e. one 
of the quasiparticle dispersions £'(T,k,q changes its sign when momentum k varies. 

Multimode FFLO-states also leave clear signatures on the noise correlations. For example, 
for the two mode FFLO-state i.e. the state where the gap is given by A(Ri) = Aq cos(q-Ri), 
the order parameter can be written as 

A(Ri) = Ao cos(q • R^) = Aq ' ^ 



k,k' 



3«q Ri 



-U 



-U 



E(cT,kq,-k-q) + e ^^-^^-^ E(cT,kq,-k+q 



Therefore the only non- vanishing expectation values in the momentum space are (c|,kCj,-k+q) 
and (c|^kCj.,-k-q)- This implies that an |-atom in the momentum state k is paired with j- 
atoms in the momentum states — k-|-q and — k — q. Because in the single mode FFLO-state 
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FIG. 4: [Color online] In this figure we show the FFLO-state noise correlation at the temper- 

ature T = 0. In (a) we show a cut in the z = plane while (b) shows the column integrated signal. 
Parameters we used were = 0.25{TT/d),qy = Qz = 0, polarization P = (n^— n|)/(n|-|-nj^) = 0.168, 
A/(2J) = 0.16, and {7/(2 J) = —1.86. In figure (a) correlations vanish in the dotted areas. We 
choose r + r' — (?it/m)q = and color-coding is again such that light colors imply high peaks and 
dark colors low. 



points r and — r + htci/m are correlated after free expansion, in the two mode FFLO- 
state r is correlated with — r + ^tq/m and — r — htci/m after free expansion. For this 
reason with the two mode FFLO-state one can see pronounced correlation peaks, when 
r + r'±htq/m = 0. Actually for the multimode FFLO-state, one could use noise correlations 
to perform a "Fourier analysis" of the periodic order parameter and probe more complicated 
spatial dependencies than those discussed here. 

In Fig. IHwe show the FFLO-state noise correlations (Gfj) at zero temperature. Fig. H] 
(a) shows the FFLO-state noise correlation in the z = plane, and Fig. H] (b) demonstrates 
the integrated FFLO noise correlation. As one can see, the Fermi surface shift creates 
asymmetry which remains also in the integrated signal. Gapless regions of the FFLO-state 
are again reflected as areas without correlation peaks and this is clearly presented in Fig. H] 
(a). 



IV. NOISE CORRELATIONS AT NON-ZERO TEMPERATURE 

Until now we have assumed zero temperature, but now we generalize our computations 
to non-zero temperatures. At non-zero temperatures it is useful to write the lowest energy 
state by using fermionic quasiparticle operators 7(T,k,q and 7].kq- Creation and annihilation 
operators can be written in terms of these operators as 

CT,k = Mk,q7T,k,q - ^k,q7{,k,q 

^^^^ 

q -k+q = Mk,q7i,k,q + ^^k,q7t k q" 
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Now the lowest energy state ansatz is given by 

I^G5) = n^U,q7kjO) (26) 

k 

and the lowest energy state is an ideal gas of quasiparticles. In other words the state in 
Eq. is the vacuum state for these quasiparticle operators. 

The non- vanishing two-operator expectation values in the densities in Eq. f ll2p at non-zero 
temperatures are 

(c{,kCT,k') = 5kk'(l«k,qr/(^T,k,q) + K'k,qP(l " /(^i,k,q))) ^^J) 
(c|,-k+qCi,-k'+q) = 5kk'(|Mk,q|V(^i,k,q) + |^^k,qP(l " /(^T,k,q))), 

where the Fermi-Dirac distribution is given by 

while non-vanishing four-operator expectation values in the noise correlations in Eq. (fT3|) 
at non-zero temperatures are given by 

(4,k4 -k'+qCi-k'+qCT,k) = |^^k,q|Vk',q|V(^T,k,q)/(^i,k',q) 

+ |Mk,q|>k',qr/(^T,k,q)(l " /(^T,k',q)) + | ^k'.q | >k,q | V (^i,k' ,q) ( 1 - /(^i,k,q)) 
+ |t^k',q|>k,qr(l - /(i5T,k'.q))(l " /(i5i,k,q)) 

+ 5kk'{kk,q|>k,qr[/(i5T,k,q)/(i5^i,k,q) + (1 " / (i^^T,k,q) ) ( 1 ' / (^^^i.k.q) )] } 

(4,k4,k'^T,k'CT,k) = |Mk,qnMk',qP/(-£^T,k,q)/(-£^T,k',q) (28) 
+ |Wk,q|>k',qr/(^T,k,q)(l " /(^i,k',q)) + | ^k'.q | >k,q | V (^T,k',q) ( 1 ' /(^i,k,q)) 
+ |^k',q|>k,q|'(l - /(i?i,k',q))(l " /(i^i.k.q)) 

+ 5kk'{kk,q|>k,qr[/(i^;T,k,q)/(i5^i,k,q) + (1 " / (i^^T,k,q) ) ( 1 " /(^^i.k.q))] 
- |«k,q|V(i?T,k,q) - l^k,qP(l - / (i^i.k.q) ) } • 

Because the densities again cancel out most of the terms, the noise correlations turn out to 
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be quite simple 



an{r,v') = A(t)'J2 



m 



m 



[/(i?T,k.q)/(i?i,kJ + (1 - /(i?T,k,q))(l - /(i?i,k,q))]|Uk,q|>k,. 



G,,{r,r') = A{tYYl 



[/(i?t,k,q)/(i?i,k,q) + (1 - /(i?T,k,q))(l " / (i^i.k,^) )] | «k,q | >k,c 




Gn(r,r') = A(t)^J] 



(29) 



[/(Et,M)/(i5i,kJ + (1 - /(i?T,k,q))(l - /(i^i,k,q))]|nk,q|>k,q| 



|Wk,q|'/(i5;T,k,q) - |t^k,q|'(l " /(i^^i.k.q)) 



Gu(r,r') = A(tr J] 



[/(Et,M)/(i5i,k,q) + (1 - /(i?T,k,q))(l - /(i^i,k,q))]|nk,q|>k,q 



+ (5(r-r')(n^(r,t)) 
, /it(k' - q) 




+ (5(r-r)(n^(r,t)). 



At non-zero temperatures the state is the BCS-state if q = and the densities are equal, 
the BP-state if q = and the densities are different. Since at finite temperatures the system 
can support polarization due to thermal effects even if both quasiparticle dispersions are 
always positive, we call the state gapless BP-state if q = and one of the quasiparticle 
dispersions changes its sign. 

Fig. [5] demonstrates why the noise correlations are a better indicator of the superfluidity 
than densities and how the BCS-state noise correlation between the components changes 
quite substantially when the temperature rises and the gap is reduced. Fig. [S] (a) shows 
a cut of the BCS-state density in the z = plane when the temperature is kBT/[2J) = 
0.196 (i.e. just below the critical temperature). Fig. O (b) shows a similar result just 
above the critical temperature. In Fig. [5] (c) and (d) we compare the column integrated 
BCS-state noise correlations between the components at the temperatures kBT/{2J) = 0.0, 
kBT/{2J) = 0.196, respectively. By comparing Fig. [5] (a) and (b), one can see that the 
densities are almost the same for the BCS-state and the normal-state. On the other hand 
in the normal-state the noise correlation is identically zero. Thermal effects are strongly 
present in the results in Fig. O (d). This can be seen from the fact that the correlations 
peak in two areas, close to the center of the figure as well as near the Fermi surface. In the 
zero temperature result in Fig. O (c) there are high peaks only near the Fermi surface. The 
reason for this difference between low and higher temperatures is that when the temperature 
rises the peak heights around the Fermi surface are reduced with reduced gap while closer 
to the center the thermal effects are relatively weaker. 

Fig. [6] demonstrates how the temperature contributes to the gapless BP-state superfiuid 
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FIG. 5: [Color online] Figure (a) shows a cut of the BCS-state density in the z = plane when the 
temperature is A;bT/(2J) = 0.196 (just below T^) while figure (b) shows a cut of the normal-state 
density in the z = plane when the temperature is kBT/{2J) = 0.20 (just above Tc). Figures 
(c) and (d) show the integrated BCS-state noise correlations between components at temperatures 
kBT/{2J) = 0.0 and kBT/{2J) = 0.196, respectively. The other parameters used were = = 
0.20, U/{2J) = -1.86, in figures (a) and (d) A/(2J) = 0.09 and in figure (c) A/(2J) = 0.35. 
Color-coding is again such that light colors imply high peaks and dark colors low. 



noise correlations. As one can see, when temperature rises the peakless areas become smaller 
and eventually vanish. The reason for this is in the temperature fluctuations, which enables 
the system to carry polarization even with gapped dispersions. These figures show also that 
the peakless areas vanish only gradually, when the temperature becomes non-zero. 

Finally, Fig. [7] demonstrates how temperature influences the noise correlations of the 
FFLO state. The states correspond to minima of the free energy at their respective tem- 
peratures. Figs. [7] (a) and (b) show an example of what happens to the correlation be- 
tween components on the 2; = plane when temperature rises from kBT/{2J) = 0.0475 to 
kBT/{2J) = 0.095. Figs. [7] (c) and (d) demonstrate the same for the column integrated 
noise correlations. As one can see from the figures, sharp areas without correlation peaks 
again disappear with increasing temperatures. However, the shift in the positions of the 
peaks persists even at non-zero temperatures. 

It might not be easy to detect the gapless states via noise correlations because at the 
temperature regime where the gapless BP-state becomes the lowest energy state, the noise 
correlations of the gapless BP-state and non-gapless BP-state appear quite similar. On the 
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FIG. 6: [Color online] This figure demonstrates how temperature contributes to the noise corre- 
lations of the gapless BP-state. In these figures the temperature gradually increases, while other 
things remain equal. The last figure corresponds to the minimum of the grand potential with polar- 
ization P = (n| — n|)/(n| +n|) = 0.40, average filling fraction (n| + nj)/2 = 0.30, A/(2J) = 0.318 
and coupling strength U /(2 J) = —3.0. We have chosen r + r' = 0. All figures show the cut in the 
z = plane. At the upper row the dotted areas are peakless. Color-coding is such that light colors 
imply high peaks and dark colors low. 



other hand, by monitoring the peak heights at 2; = (and r' = — r), one can identify the 
gapless states. The reason for this is that in the gapless BP-state the maximum peak height 
is always smaller than A{tY/8 while in the-non gapless state the maximum is always bigger 
than this. This can be seen from the fact that the term lukPl'^kP peaks on the Fermi surface 
corresponding to the average chemical potential. On the other hand, near this surface the 
peak heights are 

\u^\'\v^WiEt^)f{E^,^) + (1 - /{E^mW - /(^i,k))] 

^ ^[/(A + 5/i/2)/(A - 6f,/2) + (1 - /(A - 5^/2))(l - /(A + 5/x/2))], 
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FIG. 7: [Color online] Finite temperature noise correlations (G||) of the FFLO state. Figures 
(a) and (b) show cuts in the z = plane, while figures (c) and (d) are the column integrated 
correlations. Figures (a) and (c) were calculated at /cbT/(2J) = 0.0475 and figures (b) and (d) 
were calculated A;bT/(2J) = 0.095. In all figures we used the coupling strength U/{2J) = —1.86. 
In figures (a) and (c) P = 0.208, (n^ + n|)/2 = 0.46, = 0.26(7r/(i), qy = = 0, and A/(2J) = 
0.14. In figures (b) and (d) P = 0.24, (n| + n|)/2 = 0.37, q^ = 0.21(7r/d), qy = q^ = 0, and 
A/(2J) = 0.046. We have chosen r + r' — {ht/m)(i = 0. In figure (a) correlations vanish in the 
dotted areas. Color-coding is again such that light colors imply high peaks and dark colors low. 



where 5/i = /if — fii and then the term in the square brackets is bigger than 1/2, when 
A ± 6fi/2 > (and the state is gapped) and otherwise always lower than 1/2. 

In Fig. [8] we show the height of the maximum of peak height as a function of the polar- 
ization (Fig. M (a)) and as a function of chemical potential difference (Fig. [8] (b)). Fig. [8] (a) 
shows that when the polarization P rises then the maximum peak height becomes smaller 
and drops suddenly to zero when the gas becomes normal. With our parameters this hap- 
pens when the polarization is about 0.35. By inspecting Fig. [H] (b) it is clear that when the 
BP-state is gapless, i. e. when 6fi > 2A, the maximum peak height is indeed smaller than 
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FIG. 8: [Color online] This figure demonstrates how the polarization contributes to the maximum 
peak height, (a) shows the maximum peak height as a function of the polarization while (b) shows 
the maximum peak height as a function of the chemical potential difference — ni)/{2A) = 
Sfi/{2A). We used the parameters kBT/(2J) = 0.45 and an average filling fraction (n^ + ri|)/2 = 
0.30. We have canceled the time dependence of the maximum peak height by dividing the maximum 
peak height by the scaling factor A{t)^. The horizontal lines show value 1/8. 



1/8, a result which was made plausible above. 



V. NOISE CORRELATIONS IN A ONE-DIMENSIONAL LATTICE 



We call the lattice one- dimensional (ID) when ^ Jy = Jz- In other words a ID lattice 
can be realized with three dimensional (3D) lattice with a very strong confinement in two 
directions and weaker confinement in the third one. Since the FFLO type states are more 



favorable in ID than they are in 3D [31|, 132, 133|, 13J] and the differences between different 
phases can be more pronounced in ID than in 3D, we now consider the noise correlations 
in ID lattices. We assume that the gas in a ID lattice is released only one direction (in our 
case x-direction). The released gas creates a set of almost identical one dimensional tubes. 

Fig. [9] shows an example of the phase diagram in a ID lattice [13] • It is clear that the 
FFLO region is remarkably large and occupies the entire superfiuid region at T = 0. On the 
other hand the expected region for phase separation between the normal gas and the BCS- 
state (red) is dramatically smaller than in 3D |22i]. When the temperature rises there are 
phase transitions between the FFLO-state and, depending on polarization, phase separation, 
the BCS/BP states, or the normal gas phases. On Fig we also denote by dashed lines 
those values of temperature and polarization we have used in our following noise correlation 
computations. 

Fig. [To] demonstrates how temperature contributes to the noise correlations in a ID 
lattice. As is clear, when the temperature becomes non-zero the gapless region does not 
vanish suddenly. Also, there is a phase transition between the FFLO-state and the gapless 
BP-state when the temperature is (with our parameters) between ksT / [23^) = 0.095 and 
ksT / {2Jx) = 0.104. One can identify this phase transition from the fact that in the FFLO- 
state the noise correlation is not symmetric with respect to x = 0, while the gapless BP-state 
noise correlation is symmetric. For the two-mode FFLO-state we expect a mirror image of 
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FIG. 9: [Color online] The phase diagram of Fermi gas in an effectively ID-lattice 3j]. Colors (or 
shading): BCS/PB states=blue or dark gray, FFLO=yellow or light gray, phase separation=red 
or gray, and normal gas=white. The parameters where such that the average filling fraction 
(n-f + n|)/2 = 0.2, = O.OlEr, and U = -0.2Er. The dashed lines show the values P = 0.40 and 
kBT/{2J^) = 0.104 used in Fig. [JO] and in Fig. [H] 



the FFLO-state noise correlation shown here to appear. This mirror image would correspond 
to correlations with the pair momenta —qx- The figure which has been calculated at highest 
temperature {ksT / {2Jx) = 0.104) shows that the gapless BP-state noise correlation is again 
always below 1/8. 

In Fig. [11] we show how the polarization contributes to the noise correlations in a ID 
lattice. When the polarization is between P = 0.32 and P = 0.35 the non-gapless state 
becomes the gapless and as one can see that the maximum of the gapped BP-state noise 
correlation is always bigger than 1/8. When the polarization is between P = 0.35 to 
P = 0.48, there is a phase transition from the BP-superfiuid state to the FFLO-state which 
is visible in the pronounced asymmetry. 

In one-dimensional systems, exact solutions for the problem of interacting Fermions ex- 
ist [4^, 44, 45] . Moreover, density matrix renormalization group (DMRG) calculations can 
be applied in one dimension to provide exact numerical results. It is known that mean field 
results deviate considerably from the exact ones in one dimension. In Ref. [4^, mean field 
calculations within the one-dimensional Hubbard model (which is the case we consider) were 



compared to the exact solutions [43|, |4J, |45[. The ground state energies were found to be 
very similar, deviating at maximum a few percent. However, the mean field calculations 
were found to overestimate the excitation gap energy considerably, typically by a factor of 
1.5-3. The deviations are at their strongest at half filling. These results imply that while the 
low temperature structure in the phase diagram of figure 9 is not likely to be affected by the 
mean field approximation, the critical temperatures are probably higher than given by the 
exact solution. The dominance of the FFLO state at low temperatures and for a large range 
of polarizations, as shown by the mean field results in Figure 9, is firmly established also 



by recent exact [47|, |48| and DMRG numerical studies [35|, |49j . Note that the reference [35 



20 



kgT/(2J^)=0.0 



kgT/(2J^)=0.047 



0.25 



0.25 




0.15 



-0.25 0.25 
X [ht/(md)] 

kgT/(2J^)=0.095 



0.25 0.25 
X [ht/(md)] 

kgT/(2J^)=0.104 



0.5 




0.05 



ft<A> 



-0.5 -0.25 



X [ht/(md)] 





X [ht/(md)] 



0.25 0.5 



FIG. 10: [Color online] This figure demonstrates how temperature contributes to the noise cor- 
relations in a ID lattice. In these figures we gradually increase the temperature while keeping a 
polarization P = 0.40 and an average filling fraction (n^ + n^)/2 = 0.20 fixed. Figures in the upper 
row and bottom left show the FFLO-state noise correlations. The figure in bottom right shows the 
gapless BP-state noise correlation. In the first figure the dashed lines show the gapless region and 
in the bottom right figure the dashed line shows the value 1/8. We choose x + x' + htqx/m = in 
all these figures. In the BP-state Qx = 0. 



considers also noise correlations. It is interesting that our results for the FFLO correlations 
at zero temperature are quite similar with the more accurate computations using DMRG 
by Liischer et al. j35|. In particular, both results show an area in momentum space where 
correlation peaks are absent as well as pronounced asymmetry. In the single mode FFLO 
ansatz the signals from both ±g are naturally missing while DMRG computation is done 
for a trapped geometry and does show strong correlations corresponding to both ±g. In 
summary, our results on the structure of the noise correlations for different states in ID are 
essentially not affected by the use of mean field approximation, only the precise values of 
temperatures where these states occur are to be determined more carefully. 
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FIG. 11: [Color online] This figure shows how polarization contributes to the noise correlations in 
a ID lattice. Here we gradually increase the polarization while keeping the temperature fixed at 
ksT / {2Jx) = 0.104 and an average filling fraction at (n-f + nj)/2 = 0.20. The upper left figure 
shows the BCS-state noise correlation and upper right figure shows the gapped BP-state noise 
correlations. The bottom left figure shows the gapless BP-state noise correlation and the figure on 
the bottom right is the FFLO-state noise correlation. The dashed line shows the value 1/8. We 
choose X + x' + fitqx/m = in all these figures. In the BCS- and the BP-states qx = 0. 

VI. CONCLUSIONS 

In this paper we have presented, at the mean field level, the noise correlations of the two 
component Fermi gas in an optical lattice. We have shown that the noise correlations are a 
promising way to detect different phases in optical lattices. The different superfiuid phases 
(BCS, FFLO, and BP) can be distinguished via the noise correlations and by mapping the 
correlations more extensively, a "Fourier analysis" of the multimode FFLO-state is possible 
and could be used to reveal the structure of the more complicated periodic order parameters. 
We computed the noise correlations also at non-zero temperatures and demonstrated that 
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the differences between the correlations of different states can persist at finite temperatures 
and some quahtative features are insensitive to temperature. Also, regions of gapless quasi- 
particle dispersions can be visible in the noise correlation signals. 

Other probes also exist. By letting the gas expand freely, one can measure the momentum 
distribution by simply imaging the density distribution of the expanded cloud. From this 
momentum distribution one can, in principle, infer some important properties of the system. 
For example, one could detect the gapless BP-state in this way [50]. However, this method 



does not appear to be a promising way to detect modulated phases [35|. Spectroscopic 
means can also be considered. For example, using two-photon Bragg spectroscopy (sH . 5^ 
one can probe the systems response a t sp ecific momenta and in that way gain information 



on the quasi-particle excitations |53|, |5J, |55j. We are not aware of anyone analyzing this 



problem for fermions in optical lattices. RF-spectroscopy has been proposed [5S] for detect- 
ing a spatially nonuniform FFLO order parameter: an additional spectral peak appears due 
to quasiparticles at the nodes of the order parameter. One very appealing property of noise 
correlation signals is their sensitivity to interactions and pairing effects. If pairing does not 
take place, noise correlations between different components vanish. This is in contrast to 
other methods where the influence of the other component is less direct. Noise correlations 
can be especially useful in distinguishing between possible pseudo-gap effects and superflu- 
idity. Pseudo-gap implies non-condensed pairs, then noise correlations between momenta 
k -|- q and — k + q appear not only for a single or a few discrete values of q, but for a 
thermally distributed set of q values. The appearance of correlations only for a certain q 
signals condensation of the pairs. 

Density-density correlations of a one- dimensional system were also discussed quite ex- 
tensively and were shown to contain very clear information on the structure of the pair- 
wavefunction as well as on the quasi-particle dispersions. In one-dimensional systems, order 
parameter modulations are favored and possible problems with signals becoming smoothed 
out by column integration can be avoided. For this reason they are very attractive candidates 
to observe FFLO-type states experimentally. However, how to balance the requirements of 
sufficiently reduced dimensionality with sufficient phase coherence 33|] is to a large extent 
still unclear. 

It is physically possible that in a three-dimensional system a phase-separation occurs 



between the normal state and the BCS-state if the system is polarized [41|, |42, 157[. This 
means that the normal state exists in one part of the lattice and the BCS-state exists in 
another part of the lattice, typically in the core if a harmonic trap potential is included. 
However, other more exotic possibilities also exist 3J|. Interestingly, the phase-separation 



between the normal gas and a paired-state could be visible in the noise correlations between 
components. This follows from the fact that the noise correlations between components in 
the normal state vanish, whereas in the paired state the correlations are at their strongest 
around the Fermi momentum. In a lattice superimposed by a trap the local density, and 
therefore the local Fermi momentum, is different in different areas of the gas. This may 
allow to identify spatial phase separation and shell structures of normal and paired states 
from the freely expanded cloud where momentum has been mapped into position. 
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